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Based on the analysis of a two-orbital Hubbard model within a mean-field approach, we propose 
a mechanism for an orbital selective phase transition (OSPT) where coexistence of localized and 
itinerant electrons can be realized. We show that this OSPT exists both at and near half filling even 
in the absence of crystal field splittings or when bandwidths, orbital degeneracies and magnetic states 
are equal for both orbitals provided the orbitals have different band dispersions. Such conditions 
should be generally satisfied in many materials. We find that this OSPT is not sensitive to the 
strength of Hund's rule coupling and that heavy doping favors the collinear antiferromagnetic state 
over the OSPT. We discuss our results in relation to the iron pnictides. 
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I. INTRODUCTION 

Orbital selective phase transitions (OSPTs) leading to 
phases where localized and itinerant electrons coexist 
have attracted extensive interest from both experimen- 
talists^^ and theoreticians^^ since the observation in 
the metallic phase of Ca 2 _ :E Sr, I .Ru04 (0.2 < x < 0.5) 
of an anomalous behavior with a Curie- Weiss-like local 
spirPSl. In spite of the controversies regarding the applica- 
bility of such a proposal to real compounders', various 
mechanisms for OSPTs have been investigated, such as 
two orbitals with different bandwidth at half-filling^M2ll, 
away from half filling with crystal field splitting^, coex- 
istence of different orbital degeneracies with crystal field 
splitting at any filling 3 ^, or different magnetic states in 
different orbitals at half-filling^. 

Recently, various models based on an assumption of 
the coexistence of localized and itinerant electrons have 
been proposed in order to describe the magnetism in 
the new iron-based superconductors^HSZ! L ess wor k has 
been done on understanding the origin of such orbital se- 
lective phases (OSPs). As one possible mechanism for the 
OSPTs in the pnictides, bands with similar bandwidths 
having different intra-band Coulomb repulsion were sug- 
gested in analogy to the mechanism of a difference in 
bandwidth 38 . However, since most of the previous stud- 
ies about the origin of OSPs are focused on the paramag- 
netic (PM) state, a correct description of the magnetism 
observed in the parent compounds of most iron pnictides 
calls for a re-investigation of the mechanism responsible 
for OSPs with magnetic order. 

In this paper, a possible OSPT mechanism is proposed 
based on a simple two-dimensional (2D) two-orbital Hub- 
bard model with both orbitals having different band dis- 
persions. We solve the model in the context of mean- field 
theory (Hartree-Fock approximation (HFA)). A compar- 
ison of our results to those obtained from the dynamical 



mean-field approximation (DMFA) at and near half fill- 
ing shows that OSPT can be qualitatively captured al- 
ready at the mean-field level without taking dynamical 
fluctuations into account. The advantage of working with 
the mean-field approach is that we will be able to investi- 
gate a large variety of possible cases not easily accessible 
within the DMFA. We will show that even in the absence 
of crystal field splittings or when bandwidths, orbital de- 
generacies, magnetic states and intra-band Coulomb re- 
pulsion are equal for both orbitals, OSPTs can still occur 
at different band fillings. We will show that it is the dis- 
tinct band dispersion in both orbitals that can be identi- 
fied as the crucial ingredient for the presence of OSPTs 
with magnetic order. The mechanism we consider is in 
fact very general since usually the strength of hybridiza- 
tions between neighboring sites in different directions is 
strongly orbital-dependent in real materials, leading to 
distinct band dispersions in different orbitals. 

II. MODEL AND METHOD 

The 2D two-orbital Hubbard model is defined as 

(ij), ({*])), 1° *7 

+ (u' - 2) ^2 n^n^i -IJ ^2 S. n -S iY , (1) 

«7>7' «7>7 ; 

where tij n = t 1 (t' ) is the intra-orbital hopping integral 
between NN (NNN) sites denoted by (ij) (((ij))) with 
band indices 7 = a, (3 in units of t. U, U' and J are the 
intra-band, inter-band Coulomb interaction and Hund's 
coupling, respectively, which fulfill the rotational invari- 
ance condition U = U' + 2J. The pair-hopping term is 
ignored as it does not affect our mean-field resulta^HUl. 
c\ (ci 7<J ) creates (annihilates) an electron in orbital 7 



2 



PM 



FM 



NAF 



0.5 





FIG. 1: (Color online) Cartoons for the different magnetically 
ordered states we use in our calculations, (a) Paramagnetic 
state. The choices of sublattice and coordinate system are 
shown, (b) Ferromagnetic, (c) Neel, (d) collinear, (e) non- 
collinear, (f) bi-collinear antiferromagnetic states. 



of site i with spin a. 



is the occupation operator, 



while rii-y = riiyf + n^, and Si 7 the spin operator. 

In order to access the true ground state in a 2D sys- 
tem with hoppings up to NNN sites, the original lattice 
is divided into two sublattices A and B (see Fig. IT] (a)), 
allowing us to consider various magnetic states in uni- 
form formulation, such as the PM state, ferromagnetic 
(FM) state with momentum Qa/b.^ — (0,0) and magne- 
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and tn'JjL = m v J*&^, where 
J 



Qa/b,i = {n, it) and to' 

(BAF) state with Qa/b,i = (0,7r) andm^ - . 
non-collinear AF (NCAF) state with Qa/b.-i — ( 7T ,'k) 

= \ m l-f + im l-y\ = 

i^Efc( c L 7 t c fe+Qi7i)l with i=A or B - Tnc correspond- 
ing cartoons for different magnetic patterns are shown in 
Fig.[l](a)-(f). 



III. COMPARISONS BETWEEN THE RESULTS 
FROM HFA AND DMFA 



In order to check the validity of our mean-field calcu- 
lations, we first compare our results with those obtained 
using the DMFA 43 45 . For this comparison, the chemi- 
cal potential rather than the filling is fixed as is usually 
done in DMFA studies, and only the NAF state is al- 
lowed as required by a two-sublattice calculation within 
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FIG. 2: (Color online) Comparison of the results from 
Hartree-Fock approximation (HFA) and dynamical mean-field 
approximation (DMFA) at t a = 1, t' a = 0.6, tp = 1, t'a = 0, 
and J/U = 0.25. a) Magnetization as a function of U/t. (b)- 
(d) ((e)-(g)) show the density of states in different phases from 
DMFA (HFA). 



the DMEAP. Fig. [2] (a) shows the sublattice magnetiza- 
tion as a function of interaction U ft for the case t a = 1, 
t' a = 0.6 and tp = 1, t'g = 0. We find that while the mag- 
netic phase transition obtained from the HFA happens 
earlier than that from the DMFA and higher magneti- 
zation is detected in the HFA -indicating that dynami- 
cal fluctuations ignored in the HFA strongly suppress the 
magnetically ordered states- the variation of the magneti- 
zation with U/t obtained from the DMFA can be qualita- 
tively reproduced by the results from the HFA. Further- 
more, all the phases given from the DMFA can be quali- 
tatively captured by the HFA as shown in Figs.[2](b)-(d) 
and (e)-(g) which depict the density of states (DOS) in 
the different phases obtained from DMFA and HFA, re- 
spectively. The OSPT, where one orbital becomes local- 
ized while the other remains metallic, is clearly detected 
by both DMFA (Fig.[2](c)) and HFA (Fig.[|(f)). The re- 
sulting OSP is sandwiched between the PM metallic state 
and the NAF insulating state as seen in Fig.[2](b), (d) for 
DMFA and Fig.[2](e), (g) for HFA. The qualitative con- 
sistency between the results from HFA and DMFA imply 
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the validity of our following discussion on the OSP as well 
as on other phases in our model at the mean-field level. 
In fact, it is already known from the DMFA -where spa- 
tial fluctuations are absent- that the PM metal-insulator 
transition which is inaccessible to the HFA is precluded 
by a magnetic phase transition in the half-filled case at 
zero temperature^ which may be qualitatively described 
by the HFA. Comparing DMFA and HFA at other hop- 
ping parameters with different t' a (not shown here), we 
concluded that dynamical fluctuations play a minor role 
in the OSPT. 



IV. THE GENERAL MECHANISM FOR OSPT 

In the following we shall investigate the case of fixed 
filling at 1/2 in contrast to the case of fixed chemical po- 
tential where the filling is changed as a function of inter- 
action U /t. All the magnetically ordered states shown in 
Fig. [I] are taken into account and the ground state is the 
one with lowest total energy. Fig. [3] (a) shows the phase 
transitions happening at t a = 1, v = 0.8 and tg = 1, 
t' g = as a function of U/t. As long as U/t < 2.88, 
the ground state is a PM metal with orbital order. In a 
small interaction region of 2.88 < U/t < 3.08, an NAF 
metal with orbital order appears. Further increasing U/t 
from 3.08 up to 4.02, the a orbital exhibits NAF insulat- 
ing behavior while the /? orbital keeps the NAF metallic 
state, indicating an OSP. Orbital order disappears in this 
interaction region. At U/t > 4.02, both orbitals display 
NAF insulating behavior. 

Though we have demonstrated that the OSPT is still 
present at fixed filling in finite dimension, the mecha- 
nism for it has not yet been identified. After analyz- 
ing the noninteracting DOS, we find that several possi- 
ble mechanisms coexist, such as (i) two orbitals having 
different bandwidth with the ratio of W a /Wg — 1.3, (ii) 
the existence of orbital order due to the different band 
dispersions of the two orbitals (t' a /t a ^t' g /tp) which can 
be viewed as the existence of an effective crystal field 
splitting, and (iii) two orbitals having distinct band dis- 
persions which leads to different shapes of the noninter- 
acting partial DOS. The last effect was not considered in 
previous DMFA studies where a semicircular DOS with 
particle-hole symmetry is usually employed for all the 
orbitals. In the following, we will reveal that orbitals 
having distinct band dispersions play a crucial role in 
the OSPT with magnetic order. 

In order to figure out the essential mechanism respon- 
sible for the OSPT observed above, we will study sep- 
arately three cases, (i) We first eliminate the effect of 
different bandwidths by rescaling the hopping parame- 
ters of the a orbital from t a = 1, t' a = 0.8 to t a = 0.769, 
t' a = 0.615 so that the ratio of t'/t a = 0.8 is retained 
while the ratio of bandwidths becomes W a /Wg = 1. 
Fig.[3](b) presents the various phases as a function of U/t 
after rescaling. Though the critical points are changed 
due to the change of the total bandwidths, all the phases 



involving OSP are preserved, indicating that such an 
OSPT exists in the absence of bandwidth differences be- 
tween orbitals. 

As a second case (ii), we remove the orbital order by 
adding an effective crystal field splitting, by which the 
half- filling condition is simultaneously satisfied at U/t = 
in both orbitals. Fig. [3] (d) shows that the OSPT is 
still present in the absence of orbital order. However, 
the states with metallic behavior in both orbitals vanish 
since the Fermi level is located right at the van Hove 
singularity in the j3 orbital at U/t — 0. We have checked 
that a small t'g which shifts the van Hove singularity away 
from the Fermi level leads to the appearance of metallic 
phases in both orbitals at finite U/t. 

As a third option (iii), we eliminate the effect of or- 
bitals having distinct band dispersions but retain the dif- 
ference in bandwidth by choosing t a — 1, t' a — 0.8 and 
tg = 0.769, t'p = 0.615 which leads to f a /t a = t'^/tg = 
0.8 and W a /Wg = 1.3. As shown in Fig.[|(c), an OSP 
is precluded by NCAF states, resulting in only two suc- 
cessive phase transitions from PM metals to NCAF in- 
sulators through NCAF metals in both orbitals. Clearly, 
the OSP will be replaced by NAF insulating states in 
both orbitals at any finite U/t if we take t a — 1, t' a — 
and tg = 1.3, t'p = 0, which means a similar disper- 
sion relation t' a /t a — t'^/tg = but different bandwidth 
W a /Wg = 1.3, since the Fermi level crosses the van Hove 
singularities in both orbitals. 

Our results so far strongly point to the fact that or- 
bitals with distinct band dispersions are crucial for the 
OSPT since it is present even though all the other mech- 
anisms mentioned above are absent while different band- 
width alone will not support the existence of OSPT when 
magnetic order is considered. Fig. [3] (e) presents a phase 
diagram in the U/t-t' a plane at t a = 1, tg = 1, and 
t'p = 0. An OSP exists in a wide region of the phase di- 
agram. The phase transitions from both NAF states to 
the NCAF state and from PM metal to NAF insulator 
are of first order (solid line), otherwise second order (dot- 
ted line) . The NAF metallic state has also been detected 

in the one-band Hubbard model with NN and NNN hop- 
p m gJH_ 

V. VARIOUS EFFECTS ON THE OBSERVED 
OSPT 

Finally, we investigate various effects on the observed 
OSPT. From the phase diagram, it is obvious that we 
should discuss two cases separately: 1) t' a /t a > 1 where 
different magnetic orders like NAF and NCAF orders 
compete with each other; 2) t' a /t a < 1 where only NAF 
order occurs. In Fig. |4j we show the results at t a = 1, 
tg = 1 and t! a = 2. We first present the effect of adding 
NNN hopping t'^ . It is found that increasing t'^ favors the 
NCAF state, which squeezes the region of the NAF OSP. 
As seen in Fig.|4](b), at t' = 0.4, the region of OSP is 
smaller than that at t'^ = (Fig. [4] (a)) and at t'^ = 0.8 
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FIG. 3: (Color online) Variation of the magnetization as a function of U/t (a) at t a = 1, t' a = 0.8 and tp = 1, t'p — 0, (b) 
at t a = 0.769, t' a = 0.615 and tp = 1, t' ? = 0, (c) at t a = 1, t' a = 0.8 and ^ = 0.769, ^ = 0.615, (d) after eliminating the 
orbital order by adding an effective crystal field splitting of A = 0.798 at t a — 0.769, t' a = 0.615 and tp = 1, t'a = 0, where 
the term for crystal field splitting is written as A(n^ — n,i a ). (e) Phase diagram in U /t-t' a at t a = 1, tp = 1, t'p = 0. Here 
J/U — 0.25 and filling is 1/2. Regions of different phases are indicated by the abbreviations defined in the text. M (I) denotes 
metal (insulator). Solid and dotted lines represent first and second order phase transitions, respectively. 



the OSP completely vanishes (see Fig. |4](c)). However, 
for the case of t' a = 0.8 (not shown), the region of the 
OSP remains unchanged at t'p — 0.4, while it is reason- 
ably replaced by NCAF state at t'a = 0.8. The effect of 
Hund's rule coupling is presented in Fig. [4] (d). Com- 
pared to Fig. [I] (a) where J/U = 0.25, the region of OSP 
is enlarged at J/U — 0.0625 and a direct first-order phase 
transition from PM metals in both orbitals to the NAF 
OSP is observed instead of two successive second-order 
phase transitions through an intermediate NAF metallic 



state at J/U = 0.25. For the case of t' a = 0.8, t' /3 = 
(not shown), a similar effect of the Hund's rule coupling 
is found. 

Fig. [4] (e) shows that at a small concentration of elec- 
tronic doping of 2.5%, the OSP with NAF order is slightly 
moved to higher values of U/t and the NCAF insulating 
states existing in the undoped case is replaced by a small 
region of OSP with NCAF order which eventually become 
CAF metallic states at larger U /t. At large doping of 
20%, only two phases with PM and CAF metallic states 
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FIG. 4: (Color online) Variation of the magnetization as a function of U jt at t a = 1, t' a = 2, tp — 1. (a) t'p = 0, J/£/ = 0.25, (b) 
t'p = 0.4, J/U = 0.25, (c) t'p = 0.8, J/17 = 0.25, and (d) t'p = 0, J/U = 0.0625 at half-filling, (e) 2.5% and (f) 20% electronic 
doping at t'p = 0, J/?7 = 0.25. Regions of different phases are indicated by the abbreviations defined in the text. M (I) denotes 
metal (insulator). 



remain and the OSP vanishes as seen in Fig. [4] (f). The 
critical value of doping concentration where the OSP dis- 
appears is around 13.6%. For the case of t' a = 0.8, = 
(not shown), the OSPT also exists at 2.5% doping but is 
excluded by PM and CAF metallic states at 20% doping. 
It is interesting to notice that CAF metallic states only 
appear when the system is doped. 20% electronic doping 
is related to the filling factor in the pnictides where 6 3d 
electrons occupy 5 3d orbitals. However, after examin- 
ing various s ets of model parameters, including those for 
the pnictide d 40 * 49 *, we should emphasize that the OSPT 
disappears whenever CAF order occurs. 

VI. DISCUSSIONS AND CONCLUSIONS 

Recently, various efforts have been made in order to 
reconcile the controversies about the origin of the CAF 
phases observed in the iron pnictides. Models contain- 



ing coup led local sp ins and itinerant electrons have been 
proposec l 32 * 34 tt ^^H Experimental data has been also in- 
terprete d in t erms of a coexistence of local and itinerant 
electron ^ 51 ! 52 !. However, such a compromise doesn't seem 
to be supported by the present study of OSPT with mag- 
netic order. We find that OSPT and CAF order tend to 
avoid each other. Also, involving the inter-orbital hop- 
pings do not favor the OSP with CAF order. On the 
other hand, it is not to be expected that increasing the 
orbital degrees of freedom will dramatically change the 
situation. Furthermore, existing mechanisms proposed 
within the PM state are in conflict with the fact that the 
low temperature phases of most pnictides are magneti- 
cally ordered, and the bandwidth of different orbitals are 
almost the same. However, quantum fluctuations, espe- 
cially spatial rather than dynamical fluctuations, which 
favors paramagnetic states, may be responsible for pos- 
sible OSPT in the pnictides. 

In summary, we propose a general mechanism for 
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an OSPT in magnetically ordered states. Different or- 
bitals with different band dispersions should be quite 
widespread in real materials. Importantly, the OSPT 
according to the presented mechanism occurs in a wide 
range of model parameters, suggesting that this mecha- 
nism could be realized in nature. 
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